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"^j Abstract 

We investigate static axially symmetric black hole solutions in a four-dimensional Einstein- Yang- 
Mills-SU(2) theory with a negative cosmological constant A. These solutions approach asymptotically 
the anti-de Sitter spacetime and possess a regular event horizon. A discussion of the main properties 
O |. of the solutions and the differences with respect to the asymptotically flat case is presented. The 

mass of these configurations is computed by using a counterterm method. We note that the A = — 3 
configurations have an higher dimensional interpretation in context of d = ff supergravity. The 
^ — ■ existence of axially symmetric monopole and dyon solutions in a fixed Schwarzschild-anti-de Sitter 

background is also discussed. An exact solution of the Einstein- Yang-Mills equations is presented in 
£Sl . Appendix. 
> ■ 

oo ■ 

: 1 INTRODUCTION 

(N . 

There has been much interest in recent years in the study of black hole solutions in gravitational theories 
with a negative cosmological constant A. The initial interest in asymptotically anti-de Sitter (AAdS) 
solutions was due to the result that (sufficiently large) Schwarzschild-AdS (SAdS) black holes are thermo- 
""^5 ■ dynamically stable P^. A renewed interest in the study of AAdS black holes has appeared following the 

AdS/conformal field theory conjecture, which proposes a correspondence between physical effects associ- 
ated with gravitating fields propagating in AdS spacetime and those of a conformal field theory (CFT) on 
the boundary of AdS spacetime |21 El ■ The AAdS black hole solutions would offer the posibility of study- 
• ing the nonperturbative structure of some conformal field theories. For example, the AdSs Hawking-Page 

phase transition is interpreted as a thermal phase transition from a confining to a deconfining phase in 
the dual D — 4, N — 4 super Yang-Mills theory |3J, while the phase structure of Reissner- Nordstrom- AdS 
black holes resembles that of a van der Waals-Maxwell liquid-gas system Kl. Also, when A < 0, the 
so-called topological black holes, whose topology of the event horizon is no longer the two-sphere S may 
appear (see |Hj for reviews of the subject). 

These results motivate at least partially attempts to find new black hole solutions in AAdS spacetimes. 
Particularly interesting are solutions violating the no hair conjecture. This conjecture states that the 
only allowed characteristics of a stationary black hole are those associated with the Gauss law, such as 
mass angular momentum and U(l) charge. Black holes with hair may be useful for probing not only 
quantum gravity, but also may be interesting tests of the AdS/CFT correspondence, particularly since we 
can find such objects which are classically stable. For AAdS black holes, it has been shown recently that 
conformally coupled scalar field can be painted as hair [7|. The existence of long range Nielson-Olesen 
vortex as hair for asymptotically AdS black holes has been investigated in |SJ |U] for SAdS black hole and 
charged black strings. 

AAdS black holes solutions with SU(2) nonabelian fields have been presented in |2B[n]. The obtained 
results are strikingly different from those valid in the asymptotically flat case (nontrivial solutions exist 
for all values of A < 0). For example, regular black hole solutions exist for continuous intervals of the 
parameter space, rather than discrete points and there are configurations for which the gauge field has no 
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zeros. Regular and black hole dyon solutions without a Higgs field have also been found Insofar as 
the no-hair theorem is concerned, it has been shown that there exist stable black hole solutions in SU(2) 
Einstein- Yang-Mills (EYM) theory that are AAdS ^J[n]. However, the same system in the presence of a 
Higgs scalar field presents only configurations with very similar properties to the asymptotically-flat space 
counterparts [T2l IT31 IH] . 

The literature on AAdS solutions with nonabelian fields is continuously growing, including stability 
analyses ^| E| > the study of configurations with a NUT charge |T7] , higher dimensional counterparts JS] 
and axially symmetric generalizations ^3 |20j . Here we remark that the EYM black hole solutions with 
A < discussed in the literature are spherically symmetric or correspond to topological black holes [2"T] 
presenting the same amount of symmetry. However, for a vanishing cosmological constant, the EYM-SU(2) 
theory is known to possess also static axially symmetric finite energy black hole solutions, as exhaustively 
discussed in [221 ESj- The situation for a nonabelian field is very different from the Einstein-Maxwell theory, 
where the static black hole solution is spherically symmetric. Representing generalizations of the well 
known spherically symmetric asymptotically flat black hole solutions |24|. these static axially symmetric 
solutions are characterized by two integers, for a given event event horizon radius. These are the node 
number k of the gauge field functions and the winding number n with respect to the azimuthal angle ip 
. The static spherically symmetric solutions have winding number n — 1. Similar to the case of regular 
configurations |25j . winding numbers n > 1 leads to axially symmetric black hole solutions. Outside their 
regular event horizon, the static A = axially symmetric black hole solutions possess non-trivial magnetic 
gauge field configurations, but they carry no global magnetic charge. In a remarkable development, these 
asymptotically flat black hole solutions have been generalized to include the effects of rotation, leading to 
nonabelian counterparts of the Kerr-Newmann solution |26) . 

The regular n > 1 axially symmetric EYM configurations of Ref. [25] have been generalized for A < 
in Ref. 20 . Although some common features are present, the results in the AAdS case are rather different 
from those valid in the A = limit, in particular presenting arbitrary values for the mass and magnetic 
charge. These distinctions arise from differences that already exist in the spherically symmetric case. 

It is reasonable to suppose that the AAdS axially symmetric regular solutions discussed in Ref. [21] can 
be generalized to include a black hole event horizon inside them. The purpose of this paper is to present 
numerical arguments for the existence of this type of configurations and to analyze the properties of the 
axially symmetric EYM-SU(2) black hole solutions, in light of the existing results for the asymptotically 
flat case, discussing the points where the differences are relevant. The numerical methods used here are 
similar to the methods succesfully employed for A < axially symmetric regular solutions. 

The mass of these solutions is computed by using a counterterm method inspired by the AdS/CFT 
correspondence. Although further research is clearly necessary, at least some of these solutions, emerging 
as consistent reduction of d = 11 supergravity on a seventh dimenensional sphere |27| . may have a relevance 
in AdS/CFT context. In this case the four-dimensional cosmological constant is completely fixed by the 
gauge coupling constant and exact solutions for configurations preserving some amount of supersymmetry 
are likely to exist. 

The paper is structured as follows: in the next Section we explain the model and derive the basic 
equations, while in Section 3 we present a detailed computation of the physical quantities of the solutions 
such as mass, magnetic charge, temperature and entropy noting the relevance of a special class of solutions 
in d = 11 supergravity. In Section 4 solutions of YM equations in a fixed SAdS background are discussed. 
Some features of the axially symmetric solutions possessing a net YM electric charge in a fixed SAdS 
black hole background are also presented in this Section. The general properties of the axially symmetric 
gravitating solutions are presented in Section 5 where we show results obtained by numerical calculations. 
We give our conclusions and remarks in the final section. The Appendix contains a derivation of an exact 
solution of the EYM equations with negative cosmological constant with a planar symmetry. 
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2 GENERAL FRAMEWORK AND EQUATIONS OF MO- 
TION 



2.1 Einstein- Yang-Mills action 

We will follow in this work most of the conventions and notations used by Klcihaus and Kunz in their 
papers on asymptotically flat EYM black holes EUI2SI- 

The equations for a static, axially symmetric SU(2) gauge field coupled to Einstein gravity with a 
cosmological term have been presented in Ref. [20] • The starting point is the EYM action 

1 = J ^^16^ - 2A > - \ T ^n\ - ^ dW=/*, (1) 
where the field strength tensor is 

Fp V = d^A u - + ie [A^A V ] , (2) 

and the gauge field A^ — \r a Ap , where e and G are the Yang-Mills (YM) and the gravity coupling 
constants. The last term in (JJJ is the Hawking-Gibbons surface term [2E|, where K is the trace of the 
extrinsic curvature for the boundary dM and h is the induced metric of the boundary. Of course, this term 
does not affect the equations of motion but it is relevant for the discussion in Section 3 of the solutions' 
mass and boundary stress tensor. 

Variation of the action (JJJ with respect to the metric g^ u leads to the Einstein equations 

R v.v - + A# M „ = 8ttGT^, (3) 

where the YM stress-energy tensor is 

= 2Tr(F MQ ^ /35 Q/3 - l g ^ Fa/3 F^). (4) 

Variation with respect to the gauge field A^ leads to the YM equations 

+ie[A fl ,F» v } = 0. (5) 



2.2 Static axially symmetric ansatz and gauge condition 

As for the globally regular static solutions, we consider a line element on the form 

ds 2 = g^dafds? = -fNdt 2 + j(— + r 2 d6 2 ) + -r 2 sin 2 6d<p 2 , (6) 

where the metric functions /, m and I are only functions of r and and we note N = 1 — Ar 2 /3. Here r 
is the radial coordinate (we are interested in the region < r < oo), t is a global time coordinate, (8, ip) 
being the usual coordinates on the sphere, with < 6 < ir,0 < (p < 2ir. The metric ansatz (JSJ generalizes 
for a nonzero A the isotropic axisymmetric line element considered by Kleihaus and Kunz in [23] (a slightly 
modified version of this metric form has been considered in |Sj for asymptotically AdS black holes with 
abelian Higgs hair). The line element 10 has two Killing vectors: dt coresponding to time translation 
symmetry and d v corresponding to rotation symmetry around the z— axis (with z = rcosO, p — rsinQ). 
The event horizon resides at a surface of constant radial coordinate r = and is characterized by the 
condition /(r^) = 0, while the metric coefficients g r r,gee and g w take nonzero and finite values. 

The construction of an axially symmetric YM ansatz has been discussed by many authors starting with 
the pioneering papers by Manton 1291 and Rebbi and Rossi |30| . The most general purely magnetic axially 
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symmetric YM-SU(2) ansatz contains nine magnetic potentials and can be easily obtained in cylindrical 
coordinates x M = (p, (p, z) 

A, = \A^(p, z)r n p + l -A^\p, z)t$ + \a^{ Pi z)t?, (7) 

where the only ^-dependent terms are the SU(2) matrices (composed of the standard (n, t 2 , t 3 ) Pauli 
matrices) 

Tp = cos np t\ + srnmp t 2 , r" = — svanip t\ + cos nip t 2 , t™ = t 3 . (8) 

Transforming to spherical coordinates, it is convenient to introduce, without any loss of generality, a new 
SU(2) basis (r",r^,r™), with 

r r n = sin 6 r n p + cos r z ™ , r e n = cos 6 T n p - sin 6 r z ™ . (9) 

The general expression (J2J takes the following form in spherical coordinates 

A, = \a^{t, 6)t? + \AW{r, 9)t% + \A$\r, 0)r» (10) 

where A^/dx^ — Ardr + A^dO + A^ dp. This ansatz is axially symmetric in the sense that a rotation 
around the z— axis can be compensated by a gauge rotation C^A = with ^ being a Lie-algebra 

valued gauge function. For the ansatz ljTU|l . = ncos 6*r™/2 — n su\9tq /2. Therefore we find F w = D^W, 
where W — A v - 

Searching for axially symmetric solutions within the most general ansatz is a difficult task, and, to our 
knowledge there are no analytical or numerical results in this case. We use in this paper a reduced ansatz, 
employed also in all previous studies on EYM solutions, with five of the gauge potentials taken identically 
zero 

A(r) _ A(B) _ At) _ A6) _ _ „ 

l\ r — A r — Jig — Ji-Q — Ay — U. 

A suitable parametrization of the four nonzero components of which factorizes the trivial 0-depencence 
is 

A (r) = ^iM) j A { p = \-H 2 {r,6), A$ = -nsm6H 3 (r,6), Af = -nsin0(l - £T 4 (r, d)), (11) 

This consistent reduction of the general ansatz satisfies also some additional discrete symmetries [SU] , [22] 
(in particular the parity reflection symmetry). The explicit expression of the field strength tensor for this 
ansatz is given in [20] • To fix the residual abelian gauge invariance we choose the usual gauge condition 

[2S|-|2n] 

rd r H x - d H 2 = 0. 

The ansatz (|1U|) contains an integer n, representing the winding number with respect to the azimutal angle 
ip. While ip covers the trigonometric circle once, the fields wind n times around. Note that the usual 
spherically symmetric static ansatz corresponds to n = 1, H\ = H 3 = 0, H 2 = H4 = w(r). Similar to 
n = 1 solutions, we define the node number k by the number of nodes of the gauge fields H 2 and H4. From 
(J3J and J5J we obtain a set of seven nonlinear elliptical partial differential equations for (/, I, to, Hi) which 
should be solved numerically. These equations are presented in |2(J|. 

Within this ansatz, the energy density of the matter fields e = — T t is given by 



n- 1 



2e 2 r 4 ? 



-{rd r H 2 + cWi) 2 + ^-r-({rd r H 3 - H x H±f + {rd r H 4 + Hx{H 3 + cot6)) 2 ) (12) 
2 

~ ((deH s - 1 + cot 6H 3 + H 2 H 4 ) 2 + {d e H 4 + cot 9(H± - H 2 ) - H 2 H 3 ) 2 ^ 
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2.3 Boundary conditions 

To obtain asymptotically AdS solutions with a regular event horizon and with the proper symmetries, we 
must impose the appropriate boundary conditions, the boundaries being the horizon and radial infinity, the 
z-axis and, because of parity reflection symmetry satisfied by the matter fields, the p-axis. The boundary 
conditions at infinity and along the z— and the p-axis (i.e for 9 = 0,7r/2) agree with those of the globally 
regular solutions |2U] . 

Similar to the A = case, we impose the horizon of the black hole to reside at a surface of constant 
radial coordinate r = r/j. 

The boundary conditions satisfied by the metric functions are 

f\r=r h — m \r=r h = l\r=rh ~ (1^) 

f\r— oo — ^\r—oo ^\r— oo 1- (1^) 

We impose also the condition 

m\ g=0 = l\g= , (15) 

which is implied by the assumption of regularity on the z— axis . The boundary conditions satisfied by 
the matter functions are very similar to those satisfied in the A = case 

Hi\ r=rh = 0, d r H2\r=r h = d r Hz\r=r h = d r H4,\ r — rh = 0, (16) 

at the event horizon and 

H2\r=oo — ^4|r=oo = w 0j H-\\ r=oa = -^3^=00 = 0, (17) 

at infinity, where this time wo is a constant (with wo = ±1 for A = 0). The relations H16|) come from the 
condition F rip \ r=rh = F r g\ r=rh = imposed by the field equations, together with a special gauge choice on 
the event horizon [261 134j . For a solution with parity reflection symmetry, the boundary conditions along 
the axes are 

#l|0=O,7r/2 = H 3 \ d=0t7I /2 — 0, ,^g> 
deH 2 \e = o,Tr/2 — 9ei?4|e=o,ir/2 = 0, 

def\e=o,7r/2 = dgm\g= 0! „/2 = del\e= ,n/2 = 0. (19) 

Therefore we need to consider the solutions only in the region < 9 < ir/2. Regularity on the z— axis 
requires also 

# 2 |e=o = H4|0=o. (20) 
Dimensionless quantities are obtained by using the following rescaling 



r -> (V4^G/e)r, A -> (e 2 /47rG)A, M -> (eG/V47rG)M, (21) 
where M is the mass of the solutions (defined in the next section). 

3 PHYSICAL QUANTITIES OF THE SOLUTIONS 

3.1 Mass and boundary stress tensor 

At spatial infinity, the line element |jfj| can be written as 

ds 2 = dsl + c^dx^dx", (22) 

where are deviations from the background AdS metric ds^. Similar to the asymptotically flat case, one 
expects the value of mass to be encoded in the functions c M „. The construction of the conserved quantities 
for an asymptotically AdS spacetime was addressed for the first time in the eighties, with various aproaches 
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(see for instance Ref. The generalization of Komar's formula in this case is not straightforward 

and requires the further subtraction of a background configuration in order to render a finite result. The 
mass of the axially symmetric regular solutions discussed in |2U| has been computed in this way, by using 
the Hamiltonian formalism of Henneaux and Teitelboim |36| . 

Here we prefer to use a different approach and to follow the general procedure proposed by Balasub- 
ramanian and Kraus |37| to compute conserved quantities for a spacetime with a negative cosmological 
constant. This technique was inspired by AdS/CFT correspondence and consists in adding suitable coun- 
terterms I c t to the action. These counterterms are built up with curvature invariants of a boundary dM. 
(which is sent to infinity after the integration) and thus obviously they do not alter the bulk equations of 
motion. 

The following counterterms are sufficient to cancel divergences in four dimensions, for vacuum solutions 
with a negative cosmological constant 



let = - 



8ttG 





2 


I 


[ d 3 xV^h 




"2 R 


'dM 


I 



(23) 



Here R is the Ricci scalar for the boundary metric h. In this section we will define also I 2 = —3/ A. 

Using these counterterms one can construct a divergence-free stress tensor from the total action 
^=4uik+^surf+ict by defining 

2 SI 1 2 
T ab = -^— = —(K ab -Kh ab -jh ab + lE ab ), (24) 

where E ab is the Einstein tensor of the intrinsic metric h ab . The efficiency of this approach has been 
demonstrated in a broad range of examples, the counterterm subtraction method being developed on its 
own interest and applications. If there are matter fields on M. additional counterterms may be needed 
to regulate the action (see e.g. [HHj)- However, we find that for a SU(2) nonabelian matter content in 
four dimensions, the prescription 1)23(1 removes all divergences (a different situation is found for the five 
dimensional AAdS nonabelian solutions where the counterterm method fails and logarithmic divergences 
are presented in the total action and the expression of mass |18p. 

In deriving the explicit expression of T ab for our solutions, we make use of the asymptotic behavior of 
the metric functions 

t 1 , /i + /2sin 2 (9 1 mi + m 2 sin 2 6 1 h + 1 2 sin 2 1 

J = H 3 l-O(-r), m=l-\ - + 0{ — ), Z = H 3 1-0(— ), (25) 

where f\, ji are undetermined constants, while 

l x = mx = —, h = -^-, m 2 = —^-, (26) 

(note that we find a different power decay as compared with the A = case |251~). 
The nonvanishing components of the boundary stress tensor at large r are 

1 (Vi , 6/2., ,.\ . „{l 



t 1 fifi , 24/ 2 . 2o \, n fl 

Direct computation shows that this stress tensor is traceless. This result is expected from the AdS/CFT 
correspondence, since even dimensional bulk theories with A < are dual to odd dimensional CFTs which 
have a vanishing trace anomaly. 
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We can use this formalism to assign a mass to our EYM solutions by writing the boundary metric in 
the following form [37] 



h^dx^dx" = -N£dt 2 + a ab (dx a + N£dt){dx b + N°dt). (28) 

If £ M is a Killing vector generating an isometry of the boundary geometry, there should be an associated 
conserved charge. The conserved charge associated with time translation is the mass of spacetime 



M 



[ d 2 xy/aN s e (29) 

JO'S 



(in the absence of an electric part of the YM potential, N® = and the solutions carry no angular 
momentum). Here £ is a spacelike hypersurface from event horizon to infinity with timelike unit normal 
and e = n tl n u T flL , is the proper energy density. By using this relation, we find that the mass of our 
solutions is given by 

m = A(!A + W. (30) 



3G V 3 17 , 

The metric restricted to the boundary h ab diverges due to an infinite conformal factor r 2 // 2 . The 
background metric upon which the dual field theory resides is 

I 2 

lab = km — hab- 

r — >oo T 

For the asymptotically AdS solutions considered here, the boundary metric is 

j ab dx a dx b = l 2 (dQ 2 + sin 2 ddtp 2 ) - dt 2 , (31) 

which is just the line-element of the (2 + 1) Einstein universe. Corresponding to the boundary metric (|31|) . 
the stress-energy tensor < r a b > for the dual theory can be calculated using the following relation |3"9~] 



V=77 ab < no >= Urn V-hh ab T bc . (32) 

r — >oc 

By using this prescription, we find that the stress tensor of the field theory has the curious form 

< tu >= A diagfl, 1, —2) + B diagfl, 3, -4) sin 2 9, with A = — L(M + v^-i-), B = (33) 

where x 1 —6, x 2 — tp, x 3 — t. Here M, are continuous variables which encode the bulk parameters. 
This tensor is covariantly conserved and manifestly traceless. A winding number n > 1 of the bulk 
configurations implies /2 ^ and thus a 8— dependence of the dual theory stress tensor (although the 
boundary metric is spherically symmetric), which is a unique property of AAdS gravitating nonabelian 
configurations. This also suggests the dual theory should also contain the integer n. 

From the AdS/CFT correspondence, we expect the nonabelian hairy black holes to be described by 
some thermal states in a dual theory formulated in a (2 + 1) Einstein Universe background. The spherically 
and axially symmetric solitons will correspond to zero-temperature states in the same theory. However, 
for the EYM action Q we do not know the underlying boundary CFT. In particular, we do not know 
what the SU(2) field corresponds to in CFT language. The bulk YM fields have nothing to do, of course, 
with the nonabelian fields of the dual gauge theory. Further work in this direction will be of great interest. 

3.2 From A = -3, d = 4 EYM to d = 11 supergravity 

However, the EYM solutions with a negative cosmological constant A = —3 may have some relevance in 
AdS/CFT context (we'll work in this section with rescaled quantities as defined by l|21flh As proven in 
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|27j . for this value of the cosmological constant, an arbitrary solution {g^^A^ ) of the four dimensional 
EYM equations ©, © gives a solution of the equations of motion of the d = 11 supergravity (note also 
that the EYM system with A > can also be embedded in d — 11 supergravity |40)'l . 

Since this reduction of the eleven dimensional theory has been explained in detail in |27| . we present 
here only the basic results. The Lagrangian of the bosonic sector of d = 11 supergravity comprises the 
eleven-dimensional metric tensor §mn and a 4- index antisymmetric tensor field Fmnpq- Following |27|. 
the eleven dimensional metric ansatz reads (with fj,, v . . . denoting indices of the four dimensional space, 
M, N . . . indices of the eleventh dimensional solution, while m, n, . . . are indices on the internal space) 

ds 2 = g MN dx M dx n = g^dx^dx" + g mn (y){dy m - K {a)m A^ (x)dx»)(dy n - K {a)n A^ (x)dx l/ ), (34) 

where g^dx^dx v = e a C3) ePr\ a p is the four-dimensional metric Jfjjl, g m ndy m dy n = e p ® e q 5 pq is the metric 
on a round 7-sphere, (e a ,e p ) being the orthonormal bases on these subspaces. iO a ' m are three Killing 
vectors of the seventh dimensional internal space generating the SU(2) Lie algebra. 

The antisymmetric tensor field F can be read from (i^ )^ being a solution of the four dimensional 
YM equations) 

where e 4 = (4!)- 1 e Q/375 e Q A e* 9 A e 7 A e s , G§> = l/2e a p^F^, while = -2V m K^ a) . 

This embedding suggests the possible existence of BPS configurations with A = —3, satisfying first order 
differential equations which may lead to exact solutions. An exact solution of the field equations where 
the round two-sphere (9, ip) is replaced by a two-dimensional space of vanishing curvature is presented in 
Appendix. 



3.3 Temperature and entropy 

The zeroth law of black hole physics states that the surface gravity k is constant at the horizon of the 
black hole solutions, where 

k 2 = -(1/4)$V^)(0*0«) • (36) 

r=r h 

To evaluate k, we use the following expansions of the metric functions at the horizon 

f(r,9) = /^)(^) 2 + 0(^ X3 
m{r,9) = m2 (9)(^—^j + O ( ! ^) (37) 

Kr,9) = l 2 (9)( r -^)\o(^ 3 
\ r h J \ r h 

Since from general arguments the temperature T is proportional to the surface gravity n, T = k/(2tt), we 
obtain the relation 

T= mi_MB. (38) 

2irr h y/m 2 (0) 

Similar to the A = case, we can show, with help of the (r 9) Einstein equation which implies 

f2m 2 ,e = 2m 2 / 2 ,e, (39) 

that the temperature T, as given in l|38|l . is indeed constant. 

For the line element © , the area A of the event horizon is given by 



.1 = 2,/ d ^^fK l (40) 
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According to the usual thermodynamic arguments, the entropy S is proportional to the area A |49| 



leading to the product 



S=±, (41) 

TS= r hO-- KM [* MemOy/hW- ( 42 ) 
4G J 

Unfortunatelly, for A < we could not derive a simple Smarr-type relation, similar to that valid for 
asymptotically flat hairy black holes [221, which relates asymptotic quantities to quantities defined on 
the event horizon (this relation can also be used to test the accuracy of the numerical results). When 
integrating the Killing identity 

V a V b K a = R bc K c , (43) 
for the Killing field K a — 6" over a spacelike hypersurface E, we find the expression 



-I p desm6 r -^(Nd r f-^f) °° =A I™ dr f d$faO^^ + torG f°° dr f dBsvxf-^Tl 
2 Jo J 3 r h J rh J Q J J rh J J 

(44) 

(here we use R = 4A and = 0). To make sense on this expression we must regularize it, by extracting a 
background contribution (5^). However, this gives useful results for very simple configurations only (with 

yj—g = y—g®)- This background subtraction is of no practical use for numerical solutions like ours. 
3.4 A computation of the Euclidean action 

The expression <|4f \ for the entropy can be derived in a more rigorous way by using Euclidean quantum 
gravity arguments. Here we start by constructing the path integral [2S| 

D\g]D[9]e- a ^ (45) 

by integrating over all metrics and matter fields between some given initial and final hypersurfaces, \T> 
corresponding here to the SU(2) potentials. By analytically continuing the time coordinate t — > ir, the 
path integral formally converges, and in the leading order one obtains 

Z ~ e~ Id (46) 

where I c i is the classical action evaluated on the equations of motion of the gravity/matter system. We 
note that the considered Lorentzian solutions of the EYM equations extremize also the Euclidean action, 
t — > it having no effects at the level of the equations of motion 3 . The value of f3 is found here by 
demanding regularity of the Euclideanized manifold as r —* rh, which together with the expansion (|37fl 
and the condition (|39|l gives (3 = l/T. The physical interpretation of this formalism is that the class 
of regular stationary metrics forms an ensemble of thermodynamic systems at equilibrium temperature 
T |41| . Z has the interpretation of partition function and we can define the free energy of the system 
F = - ft' 1 logZ. Therefore 

logZ = -/3F = S-l3M, (47) 



S = /3M-I cl , (48) 

straightforwardly follows. 

3 Note that this analytical continuation becomes problematic for nonabelian solutions with an electric potential. 
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To compute I c i, we make use of the Einstein equations, replacing the R — 2 A volume term with 
2R\ — WirGTj:. For our purely magnetic ansatz, the term T\ exactly cancels the matter field lagrangean 
in the bulk action L m = —\j2Tr(F illJ F^ v \ The divergent contribution given by the surface integral term 
at infinity in R\ is also canceled by /surface + let and we arrive at the simple finite expression 



A ( 


: 2 3 A < 


8/ 2 \ 








17 J 


4G 



£ deemey/ltfj) (49) 
Replacing now in (|48|l (where M is the mass-energy computed in Section 3.1), we find 



s= igJo d0sin0 W) rl (50) 



which is one quarter of the event horizon area, as expected. 



3.5 Nonabelian charges 

Solutions of the field equations are also classified by the nonabelian electric and magnetic charges Q e and 
Qm- For purely magnetic configurations, the usual definition used by various authors (and employed also 
in this paper) is 

QM = h j ^9 Tr i pk ° T r}> ( 51 ) 

where is the dual field strength tensor. For the boundary conditions at infinity (fTTjl . we find Qm = 
n(l — ujq). However, the integral l|51|l does not intrinsically characterize a field configuration, as it is gauge 
dependent. A gauge invariant definition for the nonabelian charges has been proposed in |42| and used in 
26 for asymptotically flat nonabelian solutions 

Qe = ^ / ' dOdtp\F 8ip \, Q M = ^ j m<p\F 6v l (52) 

where the vertical bars denote the Lie- algebra norm and the integrals are evaluated as r — * 00. One can 
verify that this definition yields the absolute value of the magnetic charge as defined in l|51|) . 



4 YANG-MILLS FIELDS IN FIXED BLACK HOLE BACK- 
GROUND 

For asymptotically AdS geometries, it has been proven useful in a number of cases to consider first the 
matter system in a fixed background, as a first step toward the study of the gravitating configurations. 
These solutions preserve the basic features of the gravitating counterparts and can be obtained much more 
easily. This situation resembles the well-known case of (asymptotically-) flat monopoles and dyons, where 
the inclusion of gravity does not change the essential picture. 

In this section we present numerical arguments for the existence of nontrivial monopoles and dyons 
solutions of pure YM equations in a four dimensional SAdS black hole background, gravity being regarded 
as a fixed external field. Although being very simple, nevertheless this model appears to contain all the 
essential features of the gravitating black hole solutions. 

The background metric is given by 

dr 2 

ds 2 = — + r 2 (d0 2 + sin 2 Odip 2 ) - Fdt 2 , (53) 
F 

_ . 2M Ar 2 
r 3 
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where M is the black hole mass, and the black hole horizon Th is located at the largest root of the equation 
r/j — 2M — = 0. Note that this line element is not on the form ©; unfortunately, the SAdS 

metric does not present a simple form in "isotropic" coordinates and we should make use of the usual 
Schwarzschild coordinates. 



4.1 Monopole solutions 

4.1.1 Spherically symmetric configurations 

We start by briefly discussing the solutions obtained for n = 1 within the magnetic ansatz (flip . In this 
case H\ = H3 = 0, H2 = H4 — uj(r) and the YM equations have the simple form 



""4 



u{uj 2 -1) 2M 2Ar , 



(54) 



where the prime denotes derivative with respect to r. This problem possesses two nontrivial limits, with 
known exact solutions, noticed for the first time in The first one is found for M — > 

"= n (55) 
y/1 - A/3r 2 

describing a monopole in AdS spacetime with unit magnetic charge. Apart from this exact solution, a 
continuous family of k — finite energy configurations is found 20 I44| . Even more remarkable, the second 
one valid as A — ► describes a nonabelian hair for a Schwarzschild black hole 

(56) 



r + 3(c- 1)' 



where c = M (3 + \/3). In this case, we find an infinite sequence of excited solutions, indexed by the node 
number k (the solution 1)56(1 corresponds to k — 1), with similar properties to the well known gravitating 
counterparts |4"5)1. 

In the A < case, the numerical results show the existence of a one-parameter family of solutions 
regular at r = with the following behavior as r — > 



Ml-ArlY 



"( r ) = UJ h + z ^ r ^ T ^(r-r h ) + 0(r-r h ) 2 1 (57) 
and the asymptotic expansion at large r 



^ = c + ^+O(l), (58) 

where wq, uj\ are constants to be determined by numerical calculations and there are no restrinctions on the 
value of uiq. Since the field equations are invariant under the transformation u — > —lu, only values of LUh > 



are considered. The boundary condition l|58l) permits a non-vanishing magnetic charge Qm = 1 — u> { 

The overall picture we find for the general case is rather similar to the one described in 1111 where gravity 
is taken into account. By varying the parameter u^, a continuum of monopole solutions is obtained. The 
same general behavior is noticed for the gauge function w, with very similar branch structure. As expected, 
the properties of the solutions depends esentially on the values of A, r^. For |A| sufficiently large, there 
exist finite energy solutions for which the gauge function has no nodes. Values of ujh > 1 are always 
allowed, leading to nodeless solutions. For a fixed value of A, we obtain finite energy solutions for a finite 
number of intervals in parameter space. We noticed always the existence of at least two such interval 
around the trivial solutions u>h = 1, 0. For large enough values of |A|, the u> function remain close to its 
initial value on the event horizon, and we find only one branch of solutions. For small enough |A|, we have 
found solutions where u> crosses the r axis and a large number of branches. The arguments presented in 
[TT| for the (linear) stability of the nodeless n = 1 monopole solutions apply directly to the nongravitating 
case. 
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4.1.2 Axially symmetric configurations 

The next step is to consider the n > 1 axially symmetric generalizations of these configurations. Subject 
to the boundary conditions H13fl - I|19|) . we solve numerically for a fixed SAdS background the set of four 
YM equations implied by JSJ. 

The algorithm we have used for all axially symmetric solutions presented in this paper is only a slightly 
modified version of the approach employed by Kleihaus and Kunz in their studies of asymptotically flat 
nonabelian configurations. The field equations are first discretized on a (r, 8) grid with N r x Ng points. The 
angular coordinate 9 runs from to n/2 and the radial coordinate goes from to some large enough value 
r-max (this value is not fixed apriori, depending on the parameters r^, A; typically r max ~ 3x 10 3 -^6x lO'V/J. 
We tested for a number of solutions that the relevant quantities are insensitive to the cut off value r max . 
The grid spacing in the r— direction is non-uniform, while the values of the grid points in the angular 
direction are given by 8k — {k — l)ir/(2(Ng — 1)). Essentially, our numerical problems come from the 
region near the event horizon, where a carefull grid choice is necessary. 

In this scheme, a new radial variable is introduced which maps the semi infinite region [r^, oo) to the 
closed region [0, 1]. Our choice for this transformation was x = (r — r^)/(r + c), where c is a properly 
chosen constant in order to minimize the numerical errors. A c ^ proven to be useful for higher mass 
solutions (however in most of the computations c takes a value smaller than 20). Typical grids have sizes 
150 x 30, covering the integration region < x < 1 and < 8 < it/2. The resulting system is solved 
iteratively until convergence is achieved. All numerical calculations for axially symmetric configurations 
are performed by using the program FIDISOL (written in Fortran), based on the iterative Newton-Raphson 
method. Details of the FIDISOL code are presented in This code requests the system of nonlinear 
partial differential equations to be written in the form F(r 1 8, u, u r , ug, u r g, u rr , ugg) — 0, (where u denotes 
the unknown functions) subject to a set of boundary conditions on a rectangular domain. The user must 
deliver to FIDISOL the equations, the boundary conditions, and the Jacobian matrices for the equations 
and the boundary conditions. Starting with a guess solution, small corrections are computed until a desired 
accuracy is reached. FIDISOL automatically provides also an error estimate for each unknown function, 
which is the maximum of the discretization error divided by the maximum of the function. 

The output of the code was analysed and visualised mainly with MATHEMATICA. 

To obtain axially symmetric solutions, we start with the n — 1 solution as initial guess and increase 
the value of n slowly. The iterations converge, and repeating the procedure one obtains in this way 
solutions for arbitrary n. The physical values of n are integers. For some of the configurations, we use 
these configurations as a starting guess on a finer grid. For magnetic monopoles solutions in a fixed SAdS 
background, the typical numerical error for the functions is estimated to be lower than 10 -3 . 

The energy density of these nongravitating solutions is given by the rt-component of the energy mo- 
mentum tensor integration over all space yields their total mass-energy 



We have obtained higher winding number generalizations for every spherically symmetric configuration. 
Also, the branch structure noticed for n = 1 is retained for higher winding number solutions. These 
solutions have very similar properties with the corresponding EYM counterparts and the general picture 
we present here applies also in the Section 5. 

We have studied a large number of configurations with n — 2, 3, 4 and node number k = 0, 1, 2 
in SAdS backgrounds with = 1 and various A between 10 -3 and 10 2 . It is not easy to extract some 
general characteristic properties of the solutions, valid for every choice of the parameters (A, luo, n, k). 
However, we have found that the functions H\ and H% present always a considerable angle-dependence. 
With increasing wo, the maximal value of the H\ and H3 increase. The functions H2 and H4 present 
a small 8 dependence, although the angular dependence generally increases with Qm- The qualitative 
behavior of the gauge functions does not change by changing the value of A (the picture presented in 
Figure 10 for gravitating gauge functions applies in this case, too). 




(59) 
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As a typical example, in Figure 1 the gauge functions Hi and the energy density e are shown for 
three different solutions with node number k = 0, 1 as a function of the radial coordinate r at angles 
9 = 0, 7r/4, 7t/2. Here the winding number is n = 3 and A = —0.01, = 1; also the energy density is 
given in units 47r/e 2 . We notice that the gauge field function Hi remains nodeless and for every solution 
with Wq > 1 it takes only negative values (Hi and H 3 are zero on the axes in Fig. la and lc). The 
behavior of the energy density strongly depends on the value of Uq. For example, as seen in Figure le, for 
the k = 1 solution, the global maximum of the energy density resides on the p axis, for a r > r^. Surfaces 
of constant energy density for nongravitating solutions have a similar form to those presented in Section 
5 in the presence of gravity and we will not present them here. 

The monopole spectrum is presented in Figure 2 for a SAdS black hole with A = —3, = 1. In this 
case, the existence of one branch of solutions was noticed, and the monopole mass strongly increases with 
the magnetic charge. 



4.2 Nongravitating dyon solutions 

Following |20| , we present here numerical arguments for the existence of finite energy nongravitating YM 
dyon solutions in a fixed SAdS background. The existence of dyon solutions without a Higgs field is a 
new feature for AdS spacetime This is a consequence of the different asymptotic behavior of the AdS 
geometries as compared to the asymptotically flat case. If A > the electric part of the gauge fields is 
forbidden for static configurations ^2E3- I n order for the boundary conditions at infinity to permit the 
electric fields and maintain a finite ADM mass we have to add scalar fields to the theory. 

The YM axially symmetric ansatz ()ll|l can be generalized to include an electric part (see e.g. |26l I48j ) 

~-n n 

A t =H 5 (r,9)^+H 6 (r,9)^-. (60) 
2e 2e 

For the time translation symmetry, we choose a natural gauge such that dA/dt=0. 

For A t ^ 0, the total energy Ij59(l should be supplemented by an electric contribution 

E e = 7 / F£F ait ^g-d 3 x (61) 



4 



/oc fx/2 MM 
J drd9 sin 9 ( (rd r H 5 + ^ ) 2 + (rd r H ( 



>r h JO 

+ ^ (( 5 ^5 - H 2 H 6 ) 2 + (d e H 6 + H 2 H 5 ) 2 + n 2 (H 6 (H 3 + cot 9) + H 4 H 5 ) 2 

We can use the existence of the Killing vector d/dt which implies F it = DiA t , and the YM equations ©, 
to express the " electric mass" l|61|l as a difference of two surface integrals 0^1 

- E e = Tr( [ {D t (A t F u ^) - A t D t (F u y/=g)<Px}) = <f Tr{A t F^}dS^ - <f Tr{A t F^}dS„. (62) 

From the conditions of local finiteness of the energy density, we find that the electric potentials necessarily 
vanish on the event horizon of a static black hole (for a rotating background, the situation may be different). 
Thus, the electric mass retains only the asymptotic contribution and, similar to the regular case, a vanishing 
magnitude of the electric potentials at infinity \A t \ implies a purely magnetic solution. 

Since these axially symmetric solutions have a nonzero electric potential, we may expect them to possess 
a nonvanishing angular momentum. The total angular momentum of a nongravitating solution is given by 

J = j T^^g~d 3 x = J 2Tr{F rv F rt + F 9ip F et }^d 3 x, (63) 

and can be expressed again as a difference of two surface integrals (SUj 

J= I 2Tr{WF^}dS fl - I 2Tr{WF^ t }dS ll . (64) 

Joo Jr h 
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The relations l]6i p -l|64 p provide also an useful tests to verify the accuracy of the numerical calculations. 
The electric charge of the dyon solutions is evaluated by using the definition i|52|) . 



4.2.1 Static dyon solutions 



A possible set of boundary conditions for the electric potentials H§ , H§ is 




r=r h — 



0, H 5 \ r 



•—oo — 



^0, #6|r=QQ ^ 0, 



(65) 



and 



d9H 5 \ e=0 ^/2 — H e \g =a ^/ 2 — 



(66) 



for a solution with parity reflection symmetry. The magnetic potentials satisfy the same set of boundary 
conditions valid in purely magnetic case. 

These solutions can be considered the counterparts of the regular YM dyons in a fixed AdS geometry 
discussed in |2U|. since they satisfy the same boundary conditions at infinity and on the axes. The 
presence of a SAdS black hole will affect the properties of these solutions in the event horizon region only. 
We remark that, although locally T* ^ 0, the configurations with the symmetries implying the above 
boundary conditions along the axes have always J = 0. We use the fact that, as 9 — > it — 9 we have 
H\ — > —Hi, i?2 — * H2, -H3 — > —H3, H4 — > H4, H5 — > H5, Hq — ► — He which implies the vanishing of the 
integral Ij63(l . 

Spherically symmetric dyon solutions (discussed in |llj for the gravitating case) are found by taking 
n = 1, H$ = u(r), H2 = H4 — u)(r) and Hi — H3 = Hq = 0. The behavior of the electric potential in a 
fixed SAdS geometry is 



at large r, where uh, Uo, Ui are constants. The expansion l|57|I . H58)l for the magnetic gauge function u)(r) 
is still valid. These boundary conditions permit non- vanishing charges Qm and Qe- 

As expected, the general properties of these solutions are rather similar to the gravitating counterparts. 
There are two adjustable shooting parameters defined on the event horizon (uh, u)h)- At the horizon, u 
starts at zero and monotonically increases asymptotically to a finite value. Solutions are found for a 
continuous set of parameters and lu^; for some limiting values of these parameters, solutions blow up. 
Given (uh, 0Jh), the general behavior of the gauge functions w,U is similar to the gravitating case; there 
are also solutions with Qm — but Qe ^ 0. Similar to the soliton solutions discussed in JT], purely 
magnetic solutions are obtained by setting Uh — 0. 

When studying dyon solutions we notice the existence always of higher node (k > 1) configurations, for 
suitable values of (uh, Wfc). Also, there are solutions where uj does not cross the r axis. For a fixed value 
of u>h, the number of nodes is determined by the value of the parameter m. Typical spherically symmetric 
solutions are displayed in Fig. 3. 

The same numerical method described above is used to obtain higher winding number static dyon 
solutions in a fixed SAdS black hole background (the typical relative error is estimated to be of the order 
of 10~ 3 ). The spherically symmetric n — 1 solution gives the input data in the numerical iteration and 
it fixes the values ujq, uq in the boundary conditions at infinity. Axially symmetric n > 1 generalizations 
seem to exist for every spherically symmetric configuration. Here the inclusion of the electric potential 
does not seem to change qualitatively the general picture found in the purely magnetic case. As expected, 
for a given value of the magnetic charge, the total mass of the solutions increases with the electric charge. 
The shape of the energy density for the monopole solutions is retained for the n > 1 dyon solutions. 



u(r) = u h {r - r h ) + 0(r 3 ), 



(67) 



at the event horizon and 



u = u + — + 0( — ) 



r r 



(68) 
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As typical axially symmetrical static dyon configurations, we show in Figure 4 the radial dependence 
of the gauge functions Hi and the energy density e (in units 47r/e 2 ) for three solutions with the same 
magnetic charge and magnitude of the electric potential at infinity and n = 1, 2, 3. The results are given 
for three different angles. The value of the cosmological constant is A — —3 while the event horizon radius 
of the SAdS black hole is — 1. These solution have been obtained starting from a spherically symmetric 
configurations with the shooting parameters u>h = 0.9 and — 0.03. As seen in Fig. 4a-f, the gauge 
functions H 2 , H4, H 5 do not exhibit a strong angular dependence, while the functions Hi, H% and H e 
remain nodeless for luq > 0. The angular dependence of the matter functions increases with n and the 
location of the biggest angular spliting slightly moves further outward with n. The maximal values of the 
functions Hi, if 3 and Hq increase also with n, depending on the value of ujq. Similar results have been 
obtained for other values of A. 

4.2.2 Rotating dyon solutions 

However, the boundary conditions for the electric potentials H5 and H§ (|65[1 - <|66[1 do not exhaust all 
possibilities even for n = 1. Inspired by [T§|, we have searched for solutions with electric potentials 
satisfying a different set of boundary conditions (the magnetic potentials satisfiying the same boundary 
conditions as in the purely magnetic case) 

H 5 \ r = rh = H 6 \ r=rh =Q, H 5 \ r = 00 = Vcos6, H 6 \ r=OD = V sin 9, (69) 
deH 5 \e = o = H e \e =0 = 0, i?5|e=7r/2 = deH e \ g= i r /2 = 0, 

where V is a constant corresponding to the magnitude of the electric potential at infinity. The value 
of the electric charge can be obtained from the asymptotics of the electric potentials, since as r — > 00, 
H 5 ~ cos6*(y+(ci sin 2 # + c 2 )/r), H 6 ~ smd(V + (03 sin 2 6+c<i)/r) . We remark that, similar to the regular 
case, dyon solutions with these symmetries possess a nonzero angular momentum (since the integral (|63|l 
is nonvanishing in this case). 

As initial guess in the iteration procedure, we use the n = 1 static YM solution in fixed SAdS background 
discussed above and slowly increased the value of V . The typical relative error for the gauge functions 
is estimated to be on the order of 10~ 3 , while the relations (|62|l and (|64[1 are satisfied with a very good 
accuracy. All the solutions we present here have winding number n = 1 and have been obtained for a 
SAdS black hole with A = —3, ru = 1- However, a similar general behavior has been found for some 
other negative values of A. Several solutions with n = 2 have been also studied, possessing very similar 
properties. 

The situation here presents a number of similarities to rotating YM dyons in an AdS background 
|19j . For a given SAdS background, we found nontrivial rotating solutions for every spherically symmetric 
configuration we have considered. The solutions depend on two continous parameters: the values uiq of 
the magnetic potentials H2, H4 at infinity and the magnitude of the electric potential at infinity V. A 
nonvanishing V leads to rotating configurations, with nontrivial functions Hi, H3, H§, Hq and nonzero 
Qe- As we increase V from zero while keeping ujq fixed, a branch of solutions forms. This branch extends 
up to a maximal value of V , which depends on ujq. Along this branch, the total energy, electric charge, 
electric part of energy and the absolute value of the angular momentum increase continuously with V. 
This can be seen in Figure 5, where we present the properties of a typical branch of solutions in a SAdS 
background, for a fixed value of uiq. In this picture, the total energy E and the angular momentum J (in 
units 47r/e 2 ) as well as the electric charge Qe and the ratio E e /E are shown as a function of the parameter 
V. 

Depending on V, the energy of a rotating dyon can be several orders of magnitude greater than the 
energy of the corresponding monopole solution. We find that both E e /E and J/E 2 tend to constant values 
as V is increased. At the same time, the numerical errors start to increase, we obtain large values for both 
Qe and E, and for some V max the numerical iterations fail to converge. In this limit, the total energy and 
the electric charge diverge, while the magnetic charge takes a finite value. It is difficult to find an accurate 
value for V ma x, especially for large values of luq. Alternatively, we may keep fixed the magnitude of the 
electric potential at infinity and vary the parameter loq. 
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There are also solutions where Qm = and Qe ^ 0. A vanishing Qe implies a nonrotating, purely 
magnetic configuration. However, we find dyon solutions with vanishing total angular momentum ( J = 
for some u>o, V) which are not static (locally T* ^ 0). For the considered configurations, we have found 
that most of the angular momentum in (|64|l comes from the surface integral at infinity contribution, the 
event horizon integral taking small values. 

As a typical axially symmetrical configuration, we show in Figure 6 three-dimensional plots of the gauge 
functions Hi, the magnitude of the electric potential \A t \ — (H 2 + i?!) 1 / 2 , the local angular momentum 
T* and the energy density e for a solution with n = 1, total energy E — 6.72 (in units 47r/e 2 ), nonabelian 
charges Qm — —3.622 and Qe — 6.833 as a function of the coordinates p and z. We can see that both 
local energy and angular momentum are localized in a small region around the event horizon without being 
possible to distinguish any individual component. For the gauge function, the relevant oscillating region 
extends at least one order of magnitude beyond the event horizon. 

For some of the configurations we have considered, the surfaces of constant energy density reveal a 
curious shape, which is different from the other static solutions considered in this paper or other cases 
exhibited in the literature (see Figure 6g-i). 

In the presence of gravity, the field equations for a rotating dyon configuration are considerably more 
complicated, the metric ansatz © being supplemented with at least one extradiagonal term. However, we 
expect the gravitating counterparts of these solutions to retain the basic features discussed here. 

5 AXIALLY SYMMETRIC SOLUTIONS IN THE PRESENCE 
OF GRAVITY 

Now we come to the main subject of this paper, by including the backreaction of the YM fields on the 
black hole geometry. In this case, the equations of motion are considerably more complicated than in the 
pure YM case. This is a difficult numerical problem even for asymptotically flat geometries, where no 
exact solution is known yet. 

5.1 Spherically symmetric EYM solutions 

The spherically symmetric EYM black hole solutions presented in ^3 have been obtained for a 
Schwarzschild-like line element 

ds 2 = -^M-dt 2 + — U*= 2 + f 2 (d9 2 + sin 2 9dtp 2 ) , (70) 
p{ry Hir) 

where 

2m(f) A_ 2 

H{r) = 1 ±-L - -r 1 . 

r 3 

The spherically symmetric horizon resides at radial coordinate fh, solution of the equation f% — 2m(fft) — 
Af^/3 = 0. Since we need the spherically symmetric n = 1 solution as the starting point for the calculation 
of axially symmetric configurations, it is useful to reobtain these solutions by using a different set of 
coordinates required by the metric ansatz JJJJ. 

By imposing I — m and the metric functions / and m to be functions only of the coordinate r, the 
axially symmetric ansatz reduces to the spherically symmetric form 

ds 2 = -f(l ±r 2 )dt 2 + ^(-^- + r 2 (d8 2 + sin 2 6d^) ). (71) 

o / 1 — jf^ 

Of course, for the same configuration, the values of the event horizon radius will differ for these different 
metric parametrizations (|70|l , (|71|l . Comparation of the metric in Schwarzschild coordinates (|70() with the 
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"isotropic" metric lf7T)l yields 

/(*0 _r 2 dr 1 df /(1 _A r2) ^£W (72) 

m ( r ) r 2 ' r ./\- A r 2 y/H(F) r ' 3 p(f) 2 ' 



3 

Since the mass function m(f) is only known numerically, we have to numerically integrate the above 
relations to obtain the coordinate function r(r) and the metric functions m(r), f(r). This is a direct 
generalization of the method used in [23] for A = EYM solutions and employed also in (20j for regular 
EYM regular with negative cosmological constant. 

However, we have found more convenient to solve directly the n = 1 EYM equations by using the line 
element (|71|> . The field equations read for this parametrization 

/' 2Ar\ uj{uj 2 - 1) 



= 0, 

/"- V + + — ~^-2 K + TXr - 2^ ) -9( 2 /- 2w + T^ / + : ^)) = °'( 73 ) 
/ 2m r mr \ 27V r z J N \ 3 2m • 

3m' 4u' 2 f A / , 4m 2 \ 

1 4m — = 0, 



2m / 3 NJ r 2 N 
f' 2 : f'm' , 2/' 4/ 2 ^.,2,1 (w 2 -l) 2 \ A ( nI , 2r /J/ , m'/ 



2m r r 2 A^ \ / 

and are solved with a set of boundary conditions on the event horizon 

m = m 2 (r-r h ) 2 + 0(r-r h f, f = f 2 (r - r h f + 0(r - r h f, u = lu + 0(r - r h ) 2 . (74) 
Asymptotically we find 

m = 1 + ^ + °(^)' / = 1 + $ + °(^)' w = ^o + ^ + 0(^). (75) 

Using these boundary conditions, the equations l|73|l were integrated for = 1 and various values of 
A. Again, only values of u>h > have been considered. The value tOh = corresponds to a Reissner- 
Nordstrom-AdS solutions, while = 1 is the vacuum SAdS solution (although these solutions do not 
have a simple form in this coordinate system). 

As expected, the solutions obtained for the metric ansatz (|71|l present similar properties to those 
found for the Schwarzschild-like line-element l|70(l . Also, all the basic features of the spherically symmetric 
nongravitating solutions discussed in Section 4.1. survive in the presence of gravity. 



5.2 Properties of the axially symmetric EYM solutions 

To find axially symmetric solutions of the purely magnetic EYM equations, we have used the same nu- 
merical algorithm as for the YM solutions in a fixed SAdS background presented above. We have started 
always with a n = 1 EYM solution as initial guess and increase the value of n slowly (for a fixed ojq), 
until approaching the integer physical values. For every n, corrections are calculated successively, until the 
numerical solutions satisfies a given tolerance. The typical numerical error for the gravitating solutions is 
of the order of 1CP 3 , except for some n = 3 or k = 2 solutions and most of those with n = 4, where we 
found an error on the order of 10 -2 . This error depends on the magnetic charge and mass of the solutions. 
Axially symmetric generalizations of the n = 1 solutions with a large ratio M (wo) / 'M (ujo — 1) are difficult 
to obtain, with large numerical errors. Also, for some of the higher node solutions the convergence is 
slowed down and the errors are not small enough for the results to be reliable. A set of A = test runs 
was carried out, primarily designed to evaluate the code's ability to reproduce the Kleihaus-Kunz results. 
In this case, we have obtained a very good agreement with the results of |23j. 

The axially symmetric solutions depend on three continuous parameters (rh, A, luq) as well as two 
integers: the winding number n > 1 and the node number k. We have considered solutions with rh = 1, 
several negative values of A, the node number n = 2,3,4 and k < 2. Special care has been paid to 
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the properties of the nodeless solutions because these configurations are likely to be stable against linear 
perturbations. 

Once we have a solution, the horizon variables such as Tjj , A are calculated in a straightforward way 
from l|38(l . I|40|) . The mass of the solution is computed by using the relation l|30|) . extracting the values of 
the coefficients fi , f 2 from the asymptotics of the metric functions (we used also the asymptotic relation 
(|26|) to test the accuracy of our results). 

The behavior of the solutions is in many ways similar to that of the axially symmetric solitons discussed 
in |2U| . Again, starting from a spherically symmetric configuration we obtain higher winding number 
generalizations with many similar properties. Axially symmetric generalizations seem to exist for every 
spherically symmetric black hole solution. For a fixed winding number, the solutions can also be indexed in 
a finite number of branches classified by the mass and the non-Abelian magnetic charge. These branches 
generally follow the picture found for n = 1 (with higher values of mass, however). This is in sharp 
contrast to the A = case, where only a discrete set of solutions is found |23|. Also, the Kretschman scalar 
K = RijkiK 1 ^ 1 remains finite for every (r > r^, 9). 

It is rather difficult to find some general pattern valid for every considered configuration. Qualitatively, 
the YM field behavior is similar to that corresponding to solutions in a fixed SAdS background. We notice 
a similar shape for the functions Hi and also for the energy density. Similar to n — 1 case, as long as 
the YM energy contribution is small enough (the integral (|59|l is much smaller than the total mass), the 
YM solutions in a fixed SAdS background give a good approximation to a solution of the EYM equations. 
Again, for large enough values of |A|, the derivation of the functions H2 and if 4 from the event horizon 
value is very small. 

The gauge functions ff 2 , H 3 , H4 start always at (angle dependent) nonzero values on the event horizon. 
The general picture presented in Fig. 2 is valid in this case too. Starting with suitable n = 1 configurations, 
we find both solutions where H2, H4 cross the r axis and nodeless solutions. The gauge function H2 
is always almost spherically symmetric, while the gauge functions Hi and H3 are about one order of 
magnitude smaller than the functions H2 and H4. 

For the considered solutions, the metric functions m, /, I do not exhibit a strong angular dependence. 
These functions start with a zero value on the event horizon and approach rapidly the asymptotic values. 
We have also observed little dependence of the metric functions on the node number. The functions m 
and I have a rather similar shape, while the ratio m/l indicating the deviation from spherical symmetry is 
typically close to one, except in a region near the horizon. 

However, we found that, similar to the asymptotically flat case, the horizon is deformed. This defor- 
mation is parametrized by the ratio L e /L p , where L e is the circumference measured along the equator and 
L p is the circumference measured along the poles 



For most of the considered solutions, this derivation from spherical symmetry is very small, however, at 
the limit of the numerical errors. We have found typically L e /L p ~ 0.99 (although the situation may be 
different for other values of the event horizon radius). Similar to the static asymptotically flat case, it 
remains a challenge to construct static solutions with a large horizon deformation. 

To see the the winding number dependence for a fixed magnetic charge, we present in Fig. 7 three 
solutions with o> = 1-73 and n = 1, 2 and 3. In Figs. 7a-d the gauge field functions are shown, in Figs. 7e- 
g the metric functions, and in Fig. 7h the energy density of the matter fields. These two-dimensional plots 
exhibit the r dependence for three angles 8 = 0, ir/4 and ir/2. Note that for ujq > 1 the Hi, H3 functions 
remain nodeless (Hi and H$ are zero on the axes in Figs. 7a, c as required by the boundary conditions 
(|18|) ). As expected, the angular dependence of the matter functions increases with the winding number n 
and the location of the biggest angular splitting moves further outward. Also, we can see that for every 




(76) 



(77) 
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n, the global maximum of the energy density for these black hole solutions resides on the z— axis at the 
event horizon, while a pronounced minimum develops on the p— axis at the horizon. However, this is not 
a generic property and we expect the situation to be different for smaller values of and the same A, loq. 
In fact, even for = 1 we have found configurations with a maximum of the energy density along the p 
axis. 

In Figure 8 the mass M of this branch of black hole solutions is plotted as a function of the nonabclian 
magnetic charge Qm for various winding numbers. These solutions have been obtained for a cosmological 
constant A = —0.1. For an event horizon radius = 1, this value is of interest because the corresponding 
solutions combine basic features of the small |A| case (with a large number of branches and node number 
k > 0) and large |A| case (one branch and nodeless solutions only). For A = —0.1 and n > 1 and we find 
one branch of solutions with both nodeless (loq > 0) and one-node configurations (with loq taking negative 
values up to uj™ m ~ —0.155 and extending backwards to zero). Note also that for the physically more 
interesting value A = —3, the black hole configurations with r/j = 1 present only one branch of nodeless 
solutions, with very small variations of the gauge fields outside the event horizon. A detailed analysis of 
the static black hole solutions as well as rotating regular solutions with A = — 3, will be presented elsewhere 
in connection to the d — 11 supergravity embedding. 

Figure 9a shows a typical three-dimensional plot of the T\ component of the energy-momentum tensor 
as a function of the coordinates p and z, for a black hole solution with = 1, A = —1, Qm = 1.52, 
n = 2 and k = 0. In Figures 9b-d surfaces of constant energy density are presented for the same solution. 
As expected, the energy density is not constant at the horizon but angle dependent. Similar to the 
asymptotically flat case, the surfaces of constant energy density appear ellipsoidal for small values of e, 
being flatter at the poles than in the equatorial plane. With increasing values of e, a torus-like shape 
appears, with the horizon seen at the center of the torus. We did not observe a strong dependence of this 
picture on the value of the cosmological constant. 

To illustrate the dependence of the black hole solutions on the value of the cosmological constant, we 
show in Fig. 10 three n = 2 solutions with the same magnetic charge Qm = 1.52, obtained for different 
values of A. In Figs. lOa-d the gauge field functions are presented, in Figs. lOe-g the metric functions, 
and in Fig. lOh the energy density — T\. We can notice that, for increasing |A|, the metric functions 
approach the asymptotic values faster than in the limit of small cosmological constant, as the magnetic 
field concentrates near the event horizon. Given n, u)q, r^, by decreasing the cosmological constant from 
to — oo, the field variables decrease for fixed r. However, the angular dependence of the metric and 
matter functions does not change significantly with |A| (although the angular dependence of the matter 
variables on the event horizon increases). The peak of the energy density shifts inward with increasing |A| 
and increases in height. 

The mass-temperature diagrams for black hole monopole solutions at A = —0.1, r/i = 1 and three 
winding numbers are plotted in Figure 11. The magnetic potential at infinity ujq varies along each curve. 
The higher branches in this picture corresponds to values of loq < 1, while loq > 1 for the lower branches. 
The common point of these curves corresponds to the SAdS solution with ujq = 1 4 . Thus, it seems that 
the Hawking temperature of a EYM system appears to be suppressed relative to that of a vacuum black 
hole of equal horizon area (this result can be proven analytically for n = 1 only |2J|), while the mass 
always increases. 

We don't address here the problem of limiting solutions, which is still unclear even in the spherically 
symmetric case. Axially symmetric generalizations of the spherically symmetric configurations near the 
critical solutions are difficult to obtain, with large errors for the gauge functions Hi and H3. 

6 CONCLUSIONS 

In this paper we have presented numerical arguments that EYM theory with a negative cosmological con- 
stant possesses black hole static axially symmetric solutions. These static configurations are asymptotically 

4 For any value of rt, the solution with luq = 1 corresponds to a trivial gauge configuration with Hi = H3 = 0, H2 = Hi = 1 
i.e. Fp, y = and a SAdS geometry. 
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AdS, possesing a regular event horizon. They generalize to higher winding number the known spherically 
symmetric solutions, presenting angle-dependent fields at the horizon. 

We started by presenting arguments that YM-SU(2) theory possess solutions with nonvanishing mag- 
netic and electric charges and arbitrary winding number for a SAdS black hole fixed background. The 
spherically symmetric solutions we found have properties similar to the lower branch of their known grav- 
itating counterparts. Axially symmetric YM configurations with nonzero electric potentials have been 
considered as well, some of them presenting a nonvanishing total angular momentum. 

When including gravity, we have presented results suggesting the existence of axially symmetric EYM 
black hole solutions with negative cosmological constant. Like their regular counterparts, these black hole 
configurations have continuous values of mass and non-Abelian magnetic charge and present a branch 
structure. 

We have not considered the question of stability for higher winding number solutions. However, since 
some of the static spherically symmetric black hole solutions are stable ^U], there is all reason to believe, 
that there are also static axially symmetric black stable against small time-dependent perturbations. A 
rigorous proof is however desirable, analogous to the proof given for the spherically symmetric case. 

Axially symmetric EYM solutions with a different topology of the event horizon, generalizing for n > 1 
the known topological black holes [21] are also likely to exist. Rotating configurations should exist as 
well, leading to nonabelian counterparts of the Kerr-Newman-AdS abelian solution. Here we expect the 
situation to be more complicated as compared to the asymptotically flat case. For A ^ the boundary 
solutions in the asymptotic region are less restrictive allowing for a nonvanishing value of the electric 
potential at infinity which complicates the general picture. 

In fact we expect all known asymptotically flat configurations to generalize for a negative cosmological 
constant. It would be interesting to construct AAdS nonabelian solutions which possess only discrete 
symmetries |51j and to find the corresponding boundary stress tensor. Their interpretation in AdS/CFT 
context is a challange. 
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Appendix 

We consider spacetimes whose metric can be written locally in the form 

ds 2 = e 2C ^dr 2 + e 2A ^~ 2B ^dn 2 ~ e 2B ^dt 2 , (78) 

where dil 2 = d6 2 + f 2 (9)dip 2 is the metric on a two-dimensional surface of constant curvature 2k. The 
discrete parameter k takes the values 1, and —1 and implies the form of the function f(9) 

{sin 6*, for k = 1 
9, for fc = (79) 

sinh6*, for k = — 1. 

In these solutions, the topology of the two-dimensional space t = const., r = const, depends on the value 
of k. When k = 1, the metric takes on the familiar spherically symmetric form, for k — — 1 the (9, <p) sector 
is a space with constant negative curvature (also known as a hyperbolic plane), while for k — 0, this is a 
flat surface. For any value of k, the metric (I78|) has four Killing vectors, one timelike and three spacelike. 
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The most general expression for the appropriate SU(2) connection is obtained by using the standard 
rule for calculating the gauge potentials for any spacetime group 31 . Taking into account the symmetries 
of the line element l|78|l we find |23 

A = kg { U(r ' ^ T3dt + X(r ' ^ T3dV + (w(r ' * )Tl + ^ (r ' <)T2) M + {^f T3 + W(r ' <)T2 ~ ^ ^ Tl ) fd ^\ ' (80) 

For purely magnetic, static configurations (i.e. u = 0) it is convenient to take the \ — g au g e an d 
eliminate Q by using a residual gauge freedom. The remaining function uj depends only on the coordinate 
r. As a result, we obtain the YM curvature 

F = — (u'ndr A d6 + fuo'ndr A dip + (w 2 - k)fr 3 d9 A dw) , (81) 

where a prime denotes a derivative with respect to r. 

Inserting this ansatz into the action ifTJl. integrating and dropping the surface term, we find that the 
equations of motion can be derived from an effective action whose Lagrangian is given by (the dependence 
on the coupling constants G, e is eliminated by taking the rescaling J2U) 

L = G ^ y) ^^r-u(y), (82) 

where y l = (A,B,u>) and G ik = e 2A - B - c diag(l, -1, -2 e - 2A+2B ) and 

U = -ke B+c + Ae 2A - B+c + [to 2 - k) 2 e- 2A+3B+c . (83) 

The field equation for the variable C implies the contraint 

- e 2A ~ B - c (A' 2 - B' 2 ) + ke B+c + 2uj l2 e B ~ c - Ae 2A - B+c - (w 2 - k) 2 e~ 2A+3B+c = 0. (84) 

We remark that l|82|l allows for the reparametrization r — > f(r) which is unbroken by our ansatz. 
We are interested in finding a superpotential W such that U to satisfies the relation 

lk dW dW 
^ ^ dy l dy k ' 

which allows the first order Bogomol'nyi equations 

dr oy K 

solving also the second-order EYM equations (JSJ. 

In the absence of the YM fields, the superpotential W has the simple form 



W = y/ke 2A + e 4A - 2B , (87) 

leading to pure-AdS vacuum solutions. 

For A = —3 and vanishing curvature (k = 0) of the two-dimensional space r = const., t = const., it is 
possible to find a simple form of the superpotential 

Wo = e B w 2 +e 2A ~ B , (88) 

leading to a simple solution of the EYM equations. The line clement reads 

ds 2 = dr 2 + e {Bo e 4r ~ u 2 )dn 2 - _- J ^dt 2 , (89) 
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(where ujq, Bq are arbitrary real constants; Bq > for a solution with the right signature), while the gauge 
potential is 

w(r) = oj e- gr . (90) 

However, we find that, for every choice of the integration constants, the line element (1891 presents some 
unphysical properties. A direct computation reveals that the point r = ro (with e 4r ° = luq/Bq) is a 
curvature singularity (this can easily be seen by computing the invariant i? M ^i? A " y ), which is not hidden by 
an event horizon. As r — > ro, the metric potential gtt diverges, while the other metric functions as well as 
the gauge field remain finite. Also, the line element l|89|) approaches asymptotically the AdS background. 

Following the rules in |27j, this configuration can be uplifted to become a solution of the equations of 
motion for the d = 11 supergravity. Similar solutions are very likely to exist for k = ±1. 
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Figure la: The gauge function H\ is shown as a function of the radial coordinate r for the angles 8 = 0, 
7r/4 and tt/2 for three n = 3 axially symmetric monopole solutions in a fixed SAdS background with 
A = —0.01, r/j = 1. The magnetic charge and total mass of these YM solutions arc {Qm = —1-39, 
M = 1.09), (Q M = 3, M = 0.33) and (Q M = -11.73, M = 0.91). 
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Figure lb: Same as Fig. la for the gauge function if 2- 
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Fig. lc 
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Figure lc: Same as Fig. la for the gauge function H 3 . 
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Figure Id: Same as Fig. la for the gauge function H4. 
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Fig. le 




Figure le: Same as Fig. la for the energy density —T\ (in units 47r/e 2 ). 
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Figure 2: The total mass M (in units 47r/e 2 ) is plotted as a function of magnetic charge Qm for monopole 
solutions in a fixed Schwarzschild-anti de Sitter background with A = — 3, r/, = 1. The winding number 
n is also marked. 
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Fig. 3 




-1 ^ ' ' ' 1 

12 3 4 

logioto 

Figure 3: Typical nongravitating spherically symmetric dyon solutions in a SAdS background, for A = 
— 1, rh = 1, the same value of uuh = 0.9 and Uh = 0, 0.2, 0.4. The solution with Uh — corresponds to a 
magnetic monopole. The energy density e(r) is given in units 47r/e 2 . 
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Figure 4a: The gauge function H\ is shown as a function of the radial coordinate r at 9 = 0, 7r/4 and 
7r/2 for three different static dyon solutions. Here the solutions have winding numbers n = 1, 2, and 3, 
the same ratio Qm/k = —0.391, the same value of the electric potential at infinity u = 0.04 and masses 
M(n = 1) = 1.541, M(n = 2) = 1.714 and M(n = 3) = 2.071. 
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Fig. 4b 
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Figure 4b: Same as Fig. 4a for the 


gauge function H^- 
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Figure 4c: Same as Fig. 4a for the gauge function H 3 . 
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Fig. 4d 
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Figure 4d: Same as 


Fig. 4a for the 


gauge function H4. 
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Figure 4e: Same as Fig. 4a for the gauge function H 5 . 
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Fig. 4f 
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Figure 4f: Same as Fig. 4a for the gauge function H 6 . 
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Figure 4g: Same as Fig. 4a for the energy density e = — T\ (in units 47r/e 2 ). 
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Fig. 5 




Figure 5: The energy E and the angular momentum J (in units 47r/e 2 ) for rotating YM dyon solutions 
in fixed SAdS black hole geometry with A = — 3, = 1 are shown as a function on the parameter V for 
a fixed u>o = 2.15. Also shown are the electric charge Qe and the ratio E e /E. 

Fig. 6a 




Figure 6a: The magnetic gauge function H\, is shown as a function of the coordinates p = rs'm9, z = 
pcosO for a rotating dyon solution with too — 2.15, V = 2, n = 1. The energy and the angular momentum 
of this configuration (in units 47r/e 2 ) are E = 9.29 and J = 4.414 respectively, while Qe = 6.833. The 
results are obtained in a SAdS background with A = — 3, ru = 1- 
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Fig. 6b 




Figure 6b: Same as Fig. 6a for the gauge function H2. 
Fig. 6c 




Figure 6c: Same as Fig. 6a for the gauge function H 3 . 
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Fig. 6d 




Figure 6d: Same as Fig. 6a for the gauge function H4. 
Fig. 6e 




Figure 6e: Same as Fig. 6a for the magnitude of the electric gauge functions 
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Fig. 6f 




Figure 6f: Same as Fig. 6a for the component Tj£ of the energy momentum tensor associated with rotation. 
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Figure 6g: Same as Fig. 6a for the component T\ of the energy momentum tensor. 
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Fig. 6h 
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Figure 7a: The gauge function Hi is shown as a function of the radial coordinate r for the angles 6 = 0, 
tt/4 and tt/2. The parameters of these gravitating solutions are: n = 1, 2, and 3, k = 0, = 1, A = —0.1, 
the node number k = 0, Q M /n = -2.027, M(n = 1) = 1.541, M(n = 2) = 1.714 and M(n = 3) = 2.071. 
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Fig. 7b 
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Figure 7b: Same as Fig. 7a for the 


gauge function H%. 
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Figure 7c: Same as Fig. 7a for the gauge function H 3 . 
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Fig. 7d 
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Figure 7d: Same as Fig. 7a for the gauge function H4. 
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Figure 7e: Same as Fig. 7a for the metric function /. 
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Fig. 7f 
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Figure 7f: Same as Fig. 7a for the metric function I. 
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Figure 7g: Same as Fig. 7a for the metric function m. 
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Fig. 7f 
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Figure 7h: Same as Fig. 7a for the energy density e = — T\ of the solutions. 

Fig. 8 




Q M /n 

Figure 8: The total mass M is plotted as a function of magnetic charge Qm for black hole gravitating 
monopolc solutions at A = —0.1, r/j = 1. The winding number n is also marked. 
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Fig. 9a 




Figure 9a: The energy density e = — T' is shown as for a gravitating monopole solution with ujq — 
0.49, n — 2. Here the cosmological constant is A = — 1 and the black hole radius is = 1. 

Fig. 9b 




Figure 9b: A surface of constant energy density e = — Tjf = 0.014 are plotted for the solution of Fig. 9a. 
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Fig. 9c 




Figure 9d: Same as Fig. 9b for e = 0.03 
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Fig. 10a 
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Figure 10a: The gauge function Hi is shown as a function of the radial coordinate r for the angles 9 = 0, 
7r/4 and n/2. The results correspond to EYM black hole solutions with winding number n = 2, magnetic 
charge Qm — 1-52, horizon radius r/, = 1 and cosmological constants A = — 1, — 0.1 and —0.001. 

Fig. 10b 
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Figure 10b: Same as Fig. 10a for the gauge function H 2 - 
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Fig. 10c 
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Figure 10c: Same as Fig. 10a for the gauge function H 3 . 
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Fig. 11 
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Figure 11: The mass-temperature diagrams for black hole monopole solutions at A = —0.1, rh = 1 and 
three winding numbers. The value of the magnetic potential at infinity is the control parameter which 
varies along each curve. 
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